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1.  Introduction 


Recently  Stewart  and  Ludford  (1980)  and  Lu  and  Ludford  (1980)  have  put  for¬ 
ward  a  theory  of  one-dimensional  steady  combustion  waves  that  travel  faster  than 
justifies  the  use  of  the  combustion  approximation;  where  it  is  aBBumed  that  the 
Mach  number,  flame  speed  divided  by  a  characteristic  sound  speed,  is  vanishing¬ 
ly  small.  When  the  combustion  approximation  is  abandoned  the  momentum  equation 
is  retained  and  coupled  with  the  energy  equation.  This  coupling  is  essential 
in  the  description  of  any  combustion  phenomena  where  the  Mach  number  is  not  small. 
Among  the  topics  treated  in  these  papers  are  fast  deflagrations  and  detonations 
respectively  and  are  described  using  large  activation-energy  asymptotics.  The 
aim  is  to  lay  down  the  groundwork  for  a  theoretical  treatment  of  the  general  prob¬ 
lem  of  transition  from  deflagration  to  detonation. 

In  particular  Stewart  and  Ludford  show  that  the  solution  of  fast  deflagra¬ 
tions  has  a  very  simple  form  when  the  heat  release  during  reaction  is  small. 

This  limit  is  useful  in  the  sense  that  it  allows  explicit  formulas  to  be  devel¬ 
oped,  whereas,  otherwise  numerical  integrations  must  be  performed.  The  present 
paper  adds  to  Stewart  and  Ludford 's  original  results  and  uses  the  small  heat  re¬ 
lease  limit  of  the  results  for  large  activation  energy  to  analyze  the  combustion 
waves  discussed  in  the  previous  papers;  fast  deflagrations  and  weak  and  strong 
detonations.  An  interesting  result  from  detonation  theory,  made  explicit  in  the 
small  heat  release  limit,  is  that  the  minimum  wave  speed  of  a  detonation  is  al¬ 
most  always  greater  than  the  Chapman-Jouget  value. 

2.  The  governing  equations 

The  present  paper  will  cite  some  details  from  the  two  mentioned  papers. 
Stewart  and  Ludford  and  Lu  and  Ludford  will  be  referred  to  as  I  and  II  respec¬ 
tively.  The  equations  used  here  are  those  in  I  and  are  in  fact  valid  for  the 


description  of  both  steady  waves  of  combustion  theory,  deflagrations  and  deton¬ 
ations.  Tn  particular  p,  v,  T  and  Y  are  the  dimensionless  density,  fluid 
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velocity,  temperature  and  mass  fraction  of  the  deficient  reactant,  so  that  p, 
v,  T,  Y  -*■  1,  0,  1,  Y  ;  or  quiescent  values;  as  the  steady  wave  frame  coordin¬ 
ate  s  -*•  The  problem  of  one-dimensional  steady  combustion  has  a  special 

formulation  whenever  the  Prandlt  and  Lewis  numbers  are  set  equal  to  one  and  is 
governed  by 


yM2dV/ds  =  yM^(V-l)  +  (TV_1-l) 

(i) 

dr/ds  =  d2T/ds2  +  a\Ye~&^ 

(2) 

T  =  T  -  (y-l)M2V2/2 

0 

(3) 

T+aY=l+B+  (y-l)M^/2  with  0  =  aY 

0  — oo 

(M 

pV  =  1  . 

(5) 

In  the  present  formulation  V  and  t  are  the  independent  variables  and 

the  temperature  T  and  mass  fraction  Y  serve  as  auxiliary  variables  defined 

by  (3)  and  (U).  Equation  (5)  defines  the  density  p  and  the  pressure  has  been 

eliminated  by  use  of  the  ideal  gas  law.  The  other  parameters  that  appear,  y, 

M2,  a.  A,  0  are  explained  in  T&ble  1. 
o 

Thus  equation  (l)  and  (2)  are  to  be  solved  under  the  condition  that 

(V,t)  -*•  (1,1+(y-1)M2/2)  [corresponding  to  (T,Y)  -*■  (l,Y  )  as  s  -*■  -<*>]  and  that 
o  — 00 

the  solution  is  bounded  as  s  -*  +». 

3.  Activation-energy  asymptotics 

In  this  section  we  summarize  the  results  of  I  and  II  that  describe  the  so¬ 
lutions  for  deflagrations  and  detonations  in  the  limit  of  large  activation  energy, 
9  ->■  oo.  First  a  review  of  the  basic  properties  and  differences  of  deflagration 
and  detonation  waves  is  appropriate. 


Both  waves  represent  transition  solutions  of  equations  (l)  and  (2)  that  leave 
the  point  (V,t)  =  (1,1+(y-1)M2/2)  at  s  =  -.<*>  and  flow  into  the  point  (V,t)  * 

(V  ,t  )  at  s  =  +».  in  determining  the  fixed  points  of  equation  (l)  and  (2)  we 


find  that 


t  =  1  +  S  +  (y-l)M2/2  ,  (6) 

oo  O 

however  there  are  two  possible  choices  for  V  given  by 

00 

Vm  =  V+  =  [(YM^+1)  ±  /fl-M^)2  -  2(>+1)BM^]/(y*1)*^  •  (7) 

The  reality  of  V  requires  that  for  a  given  6  the  wave  speed  be  restricted 


so  that 


(1-M2)2 
_ o 

2(y-DSM' 


* 


Equality  refers  to  the  Chapman-Jouget  wave  speeds.  Equation  (8)  represents  a 
quadratic  in  M2  and  thus  there  are  two  possible  steady  waves 

0  <  £  M2cJ_  =  [1+(y+1)8  -  /fl+(Y-l>B]2-ll  <  1  .  (S 


corresponding  to  deflagrations  and 


1  <  K>CJ+  =  [l+(^+l)6  +  v^1+(y-1)6]2-1  11  V?  , 


corresponding  to  detonations. 

As  we  shall  see  for  deflagrations,  the  point  (V_,tw)  is  associated  with 
a  weak  deflagration  and  is  accessible,  while  (V+,tw),  corresponding  to  strong 
deflagrations  is  not.  For  detonations  (v_»Ta,)  corresponds  to  strong  detonations 
and  is  almost  always  accessible  while  (^+*T„)  corresponding  to  weak  detonations 
is  accessible  for  only  very  special  conditions. 


■ 


-  k  - 


3. a  The  asymptotic  structure  of  fast  deflagrations 


It  was  shown  in  I,  that  a  description  of  fast  deflagration  waves  for 
2  2 

0  <  M  <  M  _  was  uniformly  valid  in  the  limit  0  -*■  ®  for  Damkohler  numbers 
—  o  ocJ- 

D  of  the  form 


..2  2 

2  Mnm  2 

D  =  6  -2-j-  6  exp(6/T4)  , 

2TZ 


(11) 


where  T*  =  T^(M2)  is  a  determined  function  of  M2,  implicitly  defining  the 
2 

flame  speed  Mq  in  terms  of  the  flame  temperature  T#. 

By  expanding  all  dependent  quantities  as 


u  =  uq  +  0  1u1  +. . . 


(12) 


it  was  shown  that  the  leading  order  solutions  for  tq»  and  in  the  e_1 

expansion,  herein  denoted  by  the  zero  subscript  pould  be  constructed  by  solving 

reactionless  equations  (l)  and  (2)  subject  to  jump  conditions  across  the  flame 

sheet  located  at  s  =  0.  The  conditions  derivable  from  a  flame  sheet  analysis 

are  that  t  and  V  and  hence  dV  /ds  are  continuous  and  that  T  =  T„  at 
o  o  o  o  * 

s  =  0.  Thus 


s  2 

t  =  1  +  0e  +  (v-l)M  /2  for  s  <  0 
o  o 


and  t  =  1  +  B  +  (y-1)M  / 2  for  s  >  0 

o  o 


(13) 

(HO 


Vq  satisfies  the  differential  equations 

dV  _  ,  (y-1)M2+2  „  s 

■  v~  {Vn  ”  - 2g~)(V  -1.)  +  for  s  <  0  , 

dS  2y  Vo  °  ( y+1 )M2  °  YM2V 

o  o  o 

and 


(15) 


dV 


-rf  =  (V  -V  )(V  -V  )  for  s  >  0  . 

QS  2 V  V  O  +  O  - 

r  O 


(16) 
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The  problem  has  been  reduced  to  solving  equations  (15)  and  (l6)  under  the  condi¬ 
tion  that  V  (-»)  =  1  and  V  (0+)  =  V  (0~). 

o  o  o 

Integrating  the  latter  equation  is  straightforward  and  we  find  that  for 


s  >  0 


2Y  1 
y+1  V+-V 


-  Jin 


s  +  c  , 


(17) 


where  c  is  a  constant.  Since  V+  -  V  is  positive  Vq  ■+•  V  as  s  -*■  +“>  which 
is  appropriate  for  weak  deflagrations. 

As  noted  in  I,  we  cannot  integrate  (15)  without  further  assumption  because 

g 

the  right  hand  side  explicitly  contains  e  .  However  we  discern  some  informa¬ 
tion  by  replacing  eS  in  (15)  with  ij>.  Then  the  solution  to 


d<|>/ds  =  <j> 


(18) 


and  equation  (15)  and  (l6)  for  s  <  0  and  s  >  0  respectively  define  a  trajec¬ 
tory  In  the  (v0»$)  phase  plane.  The  trajectory  starts  at  (l,0)  at  s  » 

The  local  solutions  in  the  neighborhood  of  the  starting  point  are  easily  seen  to 
be  a  linear  combination  of 


e/[l+(y-l)Mo] 


(19) 


Thus  (1,0)  is  a  saddle  point  (since  <l)  and  we  leave  along  a  unique  inte¬ 
gral  curve.  Then  s  =  0  corresponds  to  a  point  (V#,l)  accordingly.  (See 

Fig.  1.)  Since  V  (0)=  V„(y,8,  M^)  is  determined  by  the  integration  for  s  <  0,  and 
o  w  o 

hence  as  a  function  of  the  parameters  y*8  and  Vp ,  constant  c  is  fixed  and 


T  (0)  =  T-  has  the  value 


T*  =  1  +  B  +  (y-l)M^(l-V^)/2  , 


(20) 


which  determines  the  flame  temperature  T#  as  a  unique  function  of  M  . 


si 


K 
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3.b  The  asymptotic  structure  of  detonations 

The  general  theory  of  detonations  developed  in  II  is  completely  analogous 
to  the  theory  of  fast  deflagrations  except  for  one  important  aspect;  there  is  not 
necessarily  a  unique  wave  speed  corresponding  to  an  acceptable  D.  As  before 
D  is  considered  specified  generally  as 

D  =  C02exp(e/T#)  ,  C  ~  0(1)  (21) 


where  C  and  T#  characterize  D. 

Then  the  construction  of  the  zeroth  order  solution  follows  as  before.  Equa¬ 
tions  (13)  and  (lU)  are  correct  for  t  and  V  satisfies  ( 15 )  and  (l6)  respec- 

o  o 

2 

tively.  However  for  detonations  (Mq  >  l)  the  character  of  the  initial  point 

(Vo,<j>)  =  (1,0)  is  that  of  a  source.  There  is  not  a  unique  curve,  but  in  fact 

a  one-parameter  family  of  curves  leaving  the  point.  In  particular  the  curves 

that  intersect.  <j>  =  1  at  some  Vq(o)  =  v#  can  be  assigned  a  flame  temperature 

Tn  which  serves  as  a  parametrization.  (See  Fig.  2.)  For  a  given  wave  speed, 

2  2 

_>  ^ocj+>  the  velocity  structure  in  front  of  the  flame  sheet  is  not  uniquely 
determined  as  it  is  for  deflagrations.  Instead,  for  a  given  fixed  M^,  speci¬ 
fication  of  T#  identifies  V*  uniquely  from  equation  (20)  and  hence  the  appro¬ 
priate  velocity  structure. 

The  restrictions  on  the  existence  of  detonation  structure  are  essentially 

2 

those  from  asking,  for  fixed  Mq,  what  values  of  V#  are  attainable,  and  for 
this  set  of  V#  does  there  exist  a  velocity  adjustment  behind  the  flame  (s  >  0)? 
In  II  it  was  shown  that  V#  has  the  permissible  values 


"mm 


V.  < 


Tmin 


<  V 


(22) 


■w. ■.  9  11^ 
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V*min  correST)on(is  to  the  value  of  V#  found  by  assuming  the  velocity  given 

by  the  adiabatic  shock,  which  adjusts  the  solution  from  (V  ,4>)  =  (1,0)  to 

p  p 

(V  ,0)  =  (((y-l)M  +2]/(y+l)M  ,0),  followed  by  the  integration  from  the  end  state 
s  o  o 

of  the  shock  (Vg,0)  to  The  lfiSt  integration  is  uniquely  determined 

as  (V^jO)  is  a  saddle  point  for  detonations. 

The  weak  detonation  has  =  V  and  V  *  V.  is  the  solution  for  the 
velocity  for  s>0.  In  II  it  is  shown  that  V„  =  V+  represents  the  maximum 
possible  velocity  at  the  flame  sheet  and  that  there  are  no  solutions  possible  for 


V#  >  V  .  Again  this  observation  can  be  made  simply  by  examining  the  nature  of 
the  point  (V+,l). 

Thus  for  every  M2  >  M2.. ,  V»  .  must  be  determined  by  a  numerical 

o  —  cJ+’  *mm  J 

integration  from  the  saddle  point  For  a  given  T„  then  (22)  represents 

a  complicated  restriction  on  M2  rewritten  as 


V-.jM)  <  1  - 


2[T#-(l+e)] 

(y-i)mJ 


i 


The  limit  of  small  heat  release: 


The  smallness  of  the  heat  release  may  be  due  either  to  a  small  heat  of  reac¬ 
tion  or  to  a  small  amount  of  reactant.  The  main  advantage  from  a  theoretical 
point  of  view  is  that  the  limit  allows  analytical  expressions  to  be  developed  for 
the  structure  of  deflagration  and  detonations  as  outlined  in  Section  2,  which 
otherwise  have  to  be  derived  by  numerical  integration. 

For  8  <<  1,  conditions  (9)  and  (10)  become 

0  <  M2  <  M2  T  =  1  -  /2(y+lT  61/2  +...  (2U) 

~  O  —  OCJ- 


for  deflagrations,  and 


-  8  - 


>  MocJ+  ■  1  *  ®1/2 

for  detonations,  t  is  still  given  by  (6),  where  as 

ao 

(y-l)M^+2 


(25) 


V,  = 


*  (y+iJm2 


+  •  •  •  9 


1-M 
o  o 


V  =  1  +  8/ ( 1-M  )  +. .. 

-  o 


(26) 

(27) 


for  deflagrations,  and  for  detonations  we  interchange  V+  and  V  in  the  above 
formulas.  The  apparent  non-uniformity  near  =  1  is  resolved  by  noting  that 
there  are  no  steady  solutions  possible  for 


2  2  2 

M  T  <  M  <  M  T 
ocJ-  o  ocJ+ 


(28) 


2  2 /2 

Thus  is  bounded  away  from  one  by  0(8  )  and  V+  are  uniformly  close 

to  one  as  8  0  as  required. 

Also  we  note  that  if  we  were  to  set  8=0,  the  theory  presented  so  far 

collapses  to  that  of  the  adiabatic  shock.  The  solution  of  which  is  given  by 

equation  (17)  where  V+  are  found  by  setting  8=0  appropriately.  With  the 

condition  that  V  (~°°)  =  1  there  is  a  non-constant  solution  only  if  M2  >  1. 

o  o 

In  which  case  (17)  predicts  that  V  (»)  =  [ (y-1)M2+2]/(y-1)M2.  Taylor's  classi- 

o  o  o 

cal  analysis  of  the  weak  shock  wave  assumed  that 


M  -  1  =  e  <<  1  and  V  =  1  +  ev'  +... 
o  o 


-9) 


Directly  from  equation  (8)  we  can  anticipate  Taylor-like  shock  structure  only 
when  the  heat  release  is  sufficiently  small.  We  will  see  that  this  is  the  case. 

U.a  Deflagrations 


Following  I,  for  M2  not  close  to  one,  we  write 
o 


V  =  1  +  BV'  +...  ,  T  =  1  +  (y-1)M  /2  +  8t’  , 

o  o  o 


(30) 


so  that 
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T  =  1  +  $(t'-(y-1)M  V* )  +. 
o  o 


(31) 


and 


vM2dV'/ds  =  -(l-M2)V'  +  x'  . 
o  o 


(32) 


The  solution  which  vanishes  at  s  =  -®  and  is  continuous  at  s  =  0  correspon¬ 
ding  to  deflagrations  is 


V’  =< 


[l+(y-l)M2]  ^eS  for  s  <  0  , 

[l-M2]“1{l-YM2tl+(Y-l)M^]"1e  °  ^  °) 


for  s  >  0 


(33) 


For  small  heat  release  the  flame  temperature  must  he  represented  as 


T-  =  1  +  et. 


(3U) 


and  the  condition  that  T  (0)  =  T-  then  leads  to 


t.  =  1/(1+(y-1)M~) 


(35) 


For  a  given  t„  equation  (35)  serves  to  determine  the  wave  speed,  i.e.,  M  . 

o 

However  as  we  have  indicated,  the  formula  (33)  clearly  shows  a  non-uniformity 
as  -*■  1.  From  the  expansion  of  the  Chapman-Jouget  Mach  number,  equation 
(2U)  we  see  that  to  resolve  this  difficulty  we  must  take 


,1/2 


=  1  +  oB  ,  o  <  -  /2(y+1)  <  0  . 


(36) 


We  consider  the  expansion 


1/2 

V  =  1  +  B  '  v  +  6v„  +. 
o  12 


(37) 


10  - 


and  substitution  of  (37)  into  (15)  and  (l6)  leads  to  the  conclusions  that 


v  *  0  ,  v2  =  e  /y  for  s  <  0  and 


=  0  ,  v2  =  (s+l)/y  for  s  >  0 


(38) 

(39) 


The  unboundedness  of  v2  as  s  -+  «  suggests  that  we  consider  a  change  of 
scale  for  s  >  0  (also  suggested  by  the  Taylor  shock  wave).  So  let 


61/2s  -  „ 


(Uo> 


and  then  there  is  a  nontrivial  balance  for  v^  given  by 


— i  =  lli  (v  -v  )(v  -v  ) 
dn  2y  1  +  1  - 


(Ul) 


where 


v  ±  =  -  [i  ±  X-2(y+1)/o2"i  . 


(U2) 


The  only  solution  to  equation  (Ul)  that  satisfies  the  condition  v^(0)  =  0 
is  found  to  be 


V1  = 


dn 

v  -v  e 


l-(v+/v  )e 


<5n 


where  6  =  — -^—5-  -2(y+l) 

( y+1 ) i 


(J»3) 

(hh) 


We  note  as  a  consequence  of  solution  (U3),  a  strong  deflagration  (i.e.  v^  -*■  v+ 
as  n  +0°)  is  not  a  possibility. 

Finally  for  the  Chapman -Jo uget  value  ctcj  =  -  *^2 ( y+lT >  v+  =  v_,  and  solv¬ 
ing  (Ul)  again  gives 


V1  "  ^1  {1 


) 


(1*5) 


/( y+l)/2n/y+l 


I 
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4.b  Detonations 

If  we  proceed  to  analyze  detonations  using  expansions  (30)  then  equation 
(32)  is  valid.  The  corresponding  solution  for  V’  that  vanishes  at  s  =  -*» 
and  is  continuous  at  s  =  0  is  then 


for  s  <  0  ,! 


(1*6) 


Again  the  flame  temperature  is  represented  as  in  (31*)  and 


t,  =  (l-yM^/d-Mf)  . 


(1*7) 


Thus  the  detonation  structure  described  here  has  a  unique  speed  determined  by 
the  flame  temperature  exactly  like  the  deflagration.  In  fact  the  deflagration 
and  detonation  analyzed  in  this  way  correspond  to  the  weak  deflagration  and  weak 
detonations,  i.e. 


V  -*  1  +  g/(l-M2)  +...  as  S'*00  . 

o  o 

The  fact  that  the  above  analysis  does  not  yield  a  description  of  the  strong 
detonation,  (strong  deflagration)  is  not  suprising  since  it  would  require  the 
strong  end  point  lie  close  to  V  =  1,  or 


(48) 


[(y-1)M2+21/(y+1)M^  -  1 


-  12  - 


which  is  only  true  if  M  -  1. 

Thus  to  recover  the  description  of  the  strong  detonation  it  is  again  neces¬ 
sary  to  consider  the  distinguished  limit  given  by  (36)  where 


a  _>  /2(y+l )  >  0 


(1*9) 


Expanding  V  as  in  (37)  we  find  in  particular  that 


V!  =  v#  ,  v2 


=  V»K  +  vff)s  +  T  +  V  for  s  "  0 


(50) 


where  v*  is  a  constant. 

We  note  that  for  v#  ^  0,  then  we  must  have  a  Taylor  adjustment  region  downstream 
as  well  as  the  upstream  adjustment  that  we  encountered  earlier.  We  assume  the 
stretch  (Uo)  and  the  equation  for  v^  becomes 


dv 

=  -i~  v  (v  +  —§7)  for  n  <  0 
dq  2y  1  1  y+1 

and  the  equation  (Ul)  holds  for  n  >  0. 

The  solution  of  (51)  is  given  by 


(51) 


v.  =  -fa  (i+A e-(o/>)T')'1 

1  y+1 


(52) 


Since  o  >  0,  v^  ->  0  as  q  -*■  -«>  satisfying  the  required  boundary  condition 
for  any  value  of  A  /  0.  In  fact  (52)  is  Taylor's  solution  of  the  weak  adiabatic 
shock  mentioned  at  the  end  of  section  3.  The  value  of  the  constant  A  is  not 
determined  and  its  choice  then  defines  the  strength  of  the  shock,  i.e.  the  value 
of  the  velocity  v#  at  the  flame  sheet.  Clearly  to  describe  the  shock  transi¬ 
tion  we  require 


A  >  0 


(53) 


otherwise  we  will  encounter  an  infinity  for  v^  at  some  negative  but  finite 
value  of  n.  Matching  requires  that 


v*  = 


-:'u  _  I 
y+1  1+7T  ’ 


(51*) 


* 


-  13  - 


hus  (53)  implies  that 


2c 

Y+l 


< 


<  0  . 


(55) 


The  downstream  adjustment  region,  q  >  0,  is  again  found  by  solving  (^l), 
so  that  the  solution  is  generally  given  by 

„  .  <#, 

1  l-[(v+-v1|)/(v_-v#)]e  n 

In  order  to  ensure  a  uniformly  bounded  solution  as  n  ->■  00  it  is  necessary  to  re< 
strict  v#  so  that 


(57) 


Since  v  <  v+  <  0,  equation  (57)  implies  that 


v*  < 


(58) 


We  note  that,  as  one  would  expect,  that  (56)  also  contains  the  simplest  solutions 
of  (Ul)  namely 


r,  =  v»  =  v  or  v.  for  q  >  0 

1  *  -  + 


(59) 


Finally  the  Chapman -Jouget  detonation  occurs  when 


v+  =  v_ 


“  ^  and  °CJ  *  /2^Y+iy  • 


(60) 


Its  solution  downstream,  q  >  0,  is  found  to  be 


V1  = 


-1  -  m  q 

2y  n 


>  + 


-i./E  . 

Y+1 


(61) 


-  Ik  - 


Thus  we  have  found  that  if 


(62) 


we  have  a  proper  description  of  the  velocity  structure  and  hence  the  entirety  of 
the  detonation  wave.  We  note  that  the  lower  limit  in  (62)  represents  the  maximum 
velocity  difference  that  the  adiabatic  shock  can  attain  while  the  upper  limit  repre 
sents  the  maximum  velocity  allowable  at  the  flame  sheet  such  that  there  is  a  velo¬ 
city  adjustment  region  behind  the  flame  sheet. 

Equation  (62)  then  implies  a  restriction  of  the  flame  temperature  T#  =  1  + 

1/2 

6  t#,  where  t#  is  related  to  v#  simply  by 


t«  =  -(y+1)v# 


(63) 


In  the  present  theory,  t#  is  assumed  to  be  specified  so  that  (62)  should  in  fact 
be  interpreted  as  a  restriction  on  a  (i.e.  the  Mach  number  corresponding  to  the 
wave  speed).  The  two  inequalities  expressed  in  (62)  lead  to  the  conclusion  that 
for  a  given  t* 


-  2t 


The  minimum  wave  speed  is  found  by  differentiating  the  right  hand  side  of  (6M 
and  setting  the  result  equal  to  zero.  Thus  we  find  that  when 


t#  =  (  y~1  )  /2/(  y+l )  *  cmin  ~  °cJ  =  ^(y+1)  .  (65) 

Thus  the  minimum  wave  speed  is  precisely  the  Chapman- Jouget  speed.  (See  Fig.  3.) 

Since  specification  of  T„  characterizes  the  Damkohler  number  D,  we  are 
lead  to  the  conclusion  that  in  the  context  of  the  small  heat  release  assumption, 
detonations  nearly  always  travel  faster  than  the  Chapman-Jouget  wave  speed.  And 
for  a  given  gas,  (i.e.  t#)  the  minimum  wave  speed  possible  is  nearly  always  great 

or  than  thr  Chapman-Jouget  wave  speed. 


Finally,  the  work  presented  here  can  be  extended  to  Pr,  Le  f  1.  Full  de¬ 
tails  of  this  extension  concerning  deflagrations  are  given  in  I.  We  have  carried 
out  the  work  for  deflagrations  near  the  Chapman -Jouget  wave  speed  and  for  all 
detonations  which  can  be  described  in  the  small  heat  release  limit  and  we  have  found 
the  differences  minor  enough  as  to  regard  the  discussion  presented  here  as  quite 


general . 
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TABLE  1 


Ratio  of  specific  heats. 

Mach  number  of  flame  with  respect 
to  the  undisturbed  speed. 

Nondimensional  heat  of  reaction. 

Activation  energy. 

Damkohler  number. 

_2 

a  DM 

~  Reference  mass  flux. 
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.  The  V  ,  <(>  phase  plane  shovn  for  deflagrations.  Note  that  equation  (15) 

with  dV/ds  =  0  defines  a  parabola  whose  intersection  with  $  =  0  and 

$  =  1  define  the  fixed  points  of  the  system. 

.  The  V  ,  $  phase  plane  shown  for  detonations.  Permissible  V-  lie  in 
o  ** 

a  range  <  v#  <  v+*  The  veak  detonation  goes  directly  into  V+  and 

has  no  downstream  adjustment. 

.  Shaded  region  shows  permissible  wave  speeds  o.  For  a  given  t. ,  o  .  <  o 

*  min 

where  amin  *s  nearly  always  greater  than  °cj‘ 


is4(v,^)= 
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